MRC  Technical  Summary  Report  #2114 


ON  LIMITS  OF  MULTIVARIATE  B-SPLINES. 

Wolfgang  Dahmen  and 
Charles  A.  Micchelli 


Mathematics  Research  Center 
University  of  Wisconsin-Madison 

610  Walnut  Street 
Madison.  Wisconsin  53706 

September  1980  ^ 

(Received  July  16,  1980) 


8  >9Sl 


Approved  for  public  rolooso 
Distribution  unlimited 


Sponsored  by 

0.  S.  Army  Research  Office 
p.  0.  Box  12211 
Research  Triangle  Park 
porth  Carolina  27709 


fi 


80  in  zn  056 


UNIVERSITY  OF  WISCONSIN-MADISON 
MATHEMATICS  RESEARCH  CENTER 


ON  LIMITS  OF  MULTIVARIATE  B-SPLINES 

*  ** 

Wolf  gam  Dahmen  and  Charles  A.  Micchelli 

Dedicated  to  I.  J.  Schoenberg  with  admiration  and  esteem. 

Technical  Summary  Report  #2114 
September  1980 

ABSTRACT 

In  a  definitive  series  of  papers  I.  J.  Schoenberg  with  H.  B.  Curry 
clarified  the  relationship  between  several  diverse  properties  of  distribution 
functions  on  R1 .  They  showed  that  a  distribution  function  is  a  limit  of  B- 
spline  distributions  if  and  only  if  the  reciprocal  of  its  Laplace  transform  is 
in  the  Polya-Laguerre  class.  When  the  distribution  function  corresponds  to  a 
density  A(x)  Schoenberg  showed  that  these  properties  are  equivalent  to  A 
being  a  Polya  frequency  function  or  that  the  convolution  transform  A*h  is 
variation  diminishing. 

The  purpose  of  this  paper  is  to  extend  some  of  these  properties  to  a 
multivariate  setting.  The  major  tool  in  this  investigation  is  a  notion  of 
multivariate  B-spline  which  we  have  both  studied  earlier. 
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SIGNIFICANCE  AND  EXPLANATION 


'''B-splines  occupy  a  pivotal  position  in  the  theory  of  univariate  spline 
functions.  They  offer  a  firm  bases  for  many  of  the  attractive  theoretical  and 
practical  properties  of  splines.  Therefore  any  advance  in  our  knowledge 
concerning  the  multivariate  B-spline  should  be  of  some  value. 

The  purpose  of  this  paper  is  to  investigate  another  aspect  of 
multivariate  B-splines.  We  succeed  in  identifying  those  distribution 
functions  which  are  limits  of  multivariate  B-splines  (Polya  distributions). 
This  question  was  solved  in  one  dimension  by  Curry  and  Schoenberg.  Some 
properties  of  Polya  distributions  on  Rvj  are  also  discussed.  £r 

A  surprising  consequence  of  this  work  is  the  relationship  of  the 
reciprocal  of  the  Laplace  transform  of  our  distribution  functions  and  the 
class  of  really  lineal  entire  functions  studied  more  than  25  years  ago  by 
Motzkin  and  Schoenberg. 
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ON  LIMITS  OF  MULTIVARIATE  B-SPLINES 


Wolfgang  Dahmen*  and  Charles  A.  Micchelli** 

Dedicated  to  L  J.  Schoenberg  with  admiration  and  esteem. 

Section  1.  INTRODUCTION 

In  a  definitive  series  of  papers  I.  J.  Schoenberg  [1,8,9]  with  H.  B.  Curry  [1] 
clarified  the  relationship  between  several  diverse  properties  of  distribution  functions 
on  R1 .  They  showed  that  a  distribution  function  is  a  limit  of  B-spline  distributions 
if  and  only  if  the  reciprocal  of  its  Laplace  transform  is  in  the  Polya-Laguerre  class. 
When  the  distribution  function  corresponds  to  a  density  A(x)  Schoenberg  showed  that 
these  properties  are  equivalent  to  A  being  a  Polya  frequency  function  or  that  the 
convolution  transform  A*h  is  variation  diminishing. 

The  purpose  of  this  paper  is  to  extend  some  of  these  properties  to  a  multivariate 
setting.  The  major  tool  in  this  investigation  is  a  notion  of  multivariate  B-spline 
which  we  have  both  studied  earlier  [2,6]. 

Section  2  of  this  paper  contains  some  preliminary  material  on  entire  functions  of 
affine  lineage.  For  the  most  part,  we  discovered  these  results  before  we  became  aware 
of  a  paper  by  Motzkin  and  Schoenberg  [7]  which  treats  the  same  subject.  Our  theorems 
are  slightly  stronger  than  those  stated  in  Motzkin  and  Schoenberg.  Moreover,  we  have 
patterned  our  proofs  closely  after  the  univariate  case  which  differs  from  the 
presentation  given  in  [7],  In  Section  3,  we  identify  all  distribution  functions  which 
are  limits  of  (multivariate)  B-spline  distributions.  Section  4  is  devoted  to  some  of 
the  structural  properties  of  these  density  functions.  Also,  in  this  section,  we  extend 
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some  important  identities  previously  known  only  for  the  multivariate  h-spline  and 
truncated  powers,  [2,61.  The  final  section  contains  some  results  pertaining  to  the 
analog  of  Polya  frequency  functions  and  variations  diminishing  convolution  transforms 
on  Rs. 

Section  2.  PRELIMINARIES 

In  this  section,  we  develop  the  required  background  material  concerning  Polya- 

Laguerre  functions  on  0s. 

s 

Let  Z*?  =  )  z.C.,  Z  =  (z.,...,z  ) ,  C  “  )  e  0s  denote  the  inner 

,  .  ]  1  I  S  »  8 

3*1 

product  on  0s. 

Definition  1_.  _A  function  f  ( z ) ,  Z  e  0s  _i£  the  Polya-Laguerre  clas3  Eg  provided 
that 


f(Z) 


0 

-AZ*Z+C  *Z 
e 


I  l  <zj 

3=1 


•z) 


I  I  (1  +  Cj 
v=1 


•Z)e 


-t  ■’•2 


where  £  <  ",  Z'*,^  f  R!  and  A  is  a  positive  seml-deflnlte  s  *  s  matrix 

3=1 

with  real  entries. 

We  will  usually  be  dealing  with  f  €  Eg  which  are  normalized  so  that  f(0)  -  1. 
For  this  class,  the  principal  result  which  we  will  need  is 
Theorem  1.  If  e  R8  and 


f  (z) 

n 


I  I  (1  +  i^.Z) 

3=1 


satisfies 


lim  f  ( ix)  *  f(ix),  x  e  RS 
n-w°  n 


uniformly  in  ilxll  <  a  for  some  a.  Then  f  has  an  analytic  extension  to  0s  which 


uniformly  on  compact  subsets  of  2s. 

PROOF .  Choose  any  y  e  Rs  with  lly  It  >  0.  Then  letting  a*  *  al!y»  we  have 

lim  f  (ity)  ■  f(ity) 
n+“  n 

uniformly  in  |t|  <  a' ,  t  e  R1  .  Using  Theorem  3.4,  Hirshman  and  Widder  [4],  p.  46  we 
conclude  that  the  convergence  above  holds  for  all  2  e  2  and  f(zy)  f  E1 .  Hence 


2  00  -c  (y)z 

toy)  =  e“a(y ) 2  +b( y ) z  —  (1  +  c  (y)z)e  3 
3=1  D 

where  a(y)  >  0,  Cj(y)»  b(y)  e  R1  and  Zc2(y)  <  ■».  An  easy  calculation  shows  that 


Vf  (0)*y  =  l  Cn,j*y 

n  .  , 

3=1 


and 


(Vf  (0)-y)2  -  y • Vf  (0 ) »y  -  £  (?n,j'y)2  . 
n  n  ^ 


Similarly,  we  have 


Vf (0)»y  *  b(y) 


and 

(7f(0)»y)2  -  yVf(0)y  -  Ec2(y)  +  2a(y)  . 

Since  all  the  derivatives  of  fn(zy)  at  the  origin  converge  to  the  respective 
derivatives  of  f { .?y )  we  conclude  that  lim  ^  C  ;^  =  t  exists  and  b(y)  =  £^*y. 

T”0  li 

Mote  that  the  zeros  of  fn(zy>  are  ( Cn ' 3  *y>  •  Now,  by  Hurwitz’s  theorem  the  zeros 

n.  Clearly,  the  cluster  values  of  {?n'^,y)  for 


of  f  t  zy )  are  limits  of  zeros  of  f 


fixed  y  are  on  one  hand  equal  to  {c^(y)}  and  on  the  other  hand  are  { C 3  »y }  where 
(C3}  are  the  cluster  values  of  in  Rs.  Therefore  we  may  take  c.(y)  «  C^*y 

and  so  we  see  that  0  <  a(y)  -  Ay*y  for  some  s  x  s  real  matrix  A.  Also,  since 

“  >  £c2  (y)  =  £({•*  »y)2 

for  all  y  e  Rs  we  have  Eilc^l2  <  00  which  establishes  that  f  e  Eg. 

So  far  we  have  shown  that  f  (z)  +  f(z)  pointwise  on  Rs.  Using  the  standard 

n 

normal  family  argument,  the  uniform  convergence  on  compact  subsets  ge  will  follow 
directly  from  the  uniform  boundedness  of  fn.  To  bound  fn  we  use  an  inequality  in 
[4] ,  page  44  which  yields 

n  Ip  |+4q 

If  ( z > I  <  TT  (i  +  l;n,3*zl)  <  e  n  " 

n  5-1 

where  p^  «  £cn'^*Z,  =  £(cn'^*z)2.  Clearly  both  pn,  qn  are  uniformly  bounded  on  any 

compact  set  of  £s.  Thus  fn  is  also  bounded  there  as  well  and  so  the  proof  is  complete. 

There  is  an  obvious  converse  to  Theorem  1  which  says  that  any  f  e  Eg  with 
f ( 0  >  ■  1  is  a  limit  of  polynomials  of  the  form 

P(z)  -  TT  <1  +  Cn,J*2),  Gn'j  e  RS  . 

5-1 

A  proof  of  this  fact  follows  easily  from  the  observation  that 

e!*z  =  lim  (l  +  -515)  . 

n-Ko 

The  identification  of  Eg  as  the  set  of  all  limits  of  "lineal”  polynomials  is  due  to 
Motskin  and  Schoenberg  [7]. 

Below  we  describe  a  multivariate  version  of  the  one-sided  univariate  Polya- 
Laguerre  functions.  These  functions  will  also  be  used  in  the  next  section. 

Definition  2.  We  will  say  that  DC  R1  _i£  _an  admissible  wedge  if  the  following 


conditions  holds 


a)  ty  €  W  for  all  t  >  0,  If  y  €  W 

b)  The  smallest  linear  space  of  Rs  containing  W  _is  of_  dimension  s 

c)  The  set  {zy  :  |z|  <  1,  !  f  f,  y  f  »l  is  a  uniqueness  set  for  entire  functions 


on  jfs. 


The  polar  set  of  W  is  defined  as 


0  s 

w  »  {x  i  x  e  R  ,  y*x  >  0,  for  all  yew}. 

Example  1 .  R®  “  {y  :  y  =  (y,j .  •  • .  ,ys> ,  y^  >  0}  is  an  admissible  wedge  which  is  self- 
polar. 

Definition  3.  We  will  call  P  a  W  -  lineal  polynomial  if  it  has  the  form 


P<  Z)  -  TT  Sj‘Z 
j-1 


where  e  w  .  ftny  function  which  is  a  limit  (uniform  on  every  compact  subset  of  £s) 


of  w  -  lineal 


is  called  a  W  -  the  lineal  entire  function. 


Our  intention  is  to  show  that  f  is  W  -  lineal  if  and  only  if  it  has  the  form 


7  T  m  -1  . 

f ( z )  -  e  5  I  I  <z3*z)  I  I  (1  +  C 3*z) 

J-i  3-1 

where  .  Z^  are  in  W°  and  J  <  »>.  Let  us  denote  this  class  by  Eg(W). 

n  j  -in 

Theorem  2.  Let  f  (z)  *  II  (1  +  Cn,3,z).  c"'3  e  W°  and  suppose  w  is  an  admissible 

", 

wedge.  If 


lim  f  (ix)  ■  f(lx),  x  f  R 
n-**»  n 


uniformly  in  Ixl  <  a  for  some  a.  Then  f  has  an  analytic  extension  to  £s  which 

_is  E#(w)  and  fn(z)  converges  uniformly  to  f  ( ;>  on  compact  subsets  of  (Ts. 

PROOF.  Choose  any  y  *  W  with  lyl  >  0  then  lim  f  (ity)  -  f(ity)  uniformly  in 

n** 

|tl  <  a'  ■  a/lyl.  Since  fn<ty)  is  a  polynomial  with  only  negative  zeros  we  know 
that  (c.f.  Karlin  (5] )  that 


-5- 


lim  f  (zy)  =  f(zy) 
n 

TV*» 


uniformly  on  compact  subsets  of  £  and  f(zy)  has  the  form 


f  ( zy  1  =  ea(y)Z  11(1+  b.(y)z) 
3=1  D 


for  some  constants  satisfying  a(y)  >  0,  b_.(y)  >  0  and  Eb^(y)  < 

Just  as  in  the  proof  of  Theorem  1,  we  may  choose  b^(y)  =  C^*y  where  tp  are  the 


cluster  values  of  {i;11'3}.  Moreover,  since  ?J*y  >  0  for  all  yew  we  have 
p  e  w°.  Also, 


^  *y  >  0  fc 


(1) 


[  Ic^’yl  <  " 


for  yew  implies  that  £  It^'yl  <  "  for  any  y  in  the  linear  span  of  w.  Thus  (1) 

1 

holds  for  all  y  e  Rs  and  it  follows  that 


l  II ;  ^  II  <  *  . 


Also,  since 


lim  Vf  (0)*y  -  7f(0)»y 
n-*»  n 


we  obtain  for  every  yew 


lim  l  'v  *  l  Cj,y  +  a(y>  • 

n-x»  j=1  j=1 


Since  W  is  an  admissible  this  equation  extends  to  all  yep5  and  no  we  conclude 

that  aty)  =  ;^*y  for  some  5*^  e  ,  £  IC^U  <  and 

3 


f(z)  =  eZ*?  II  (1  +  rP'z)  . 
1=1 


-e- 


It  remains  to  show  that  f  *  f  uniformly  on  compact  subsets  of  £s.  In  view  of 

n 

property  c)  in  Definition  2  it  is  sufficient  to  prove  f  is  bounded  on  compact 
subsets  of  If8.  To  this  end,  observe  that  for  any  yew  there  is  a  constant  M(y) 
such  that 


l  ! :n'^ *y |  *  J  ;n'^*y  <  M(y) 

1  1 


for  all  n.  Since  W  is  an  admissible  wedge  it  follows  that 


l  lljn'^lt  <  M 

j“1 


for  some  constant  M  and  consequently 

n 

l  icn°ii«z« 

If  (Z)|  »  ITT  (1  +  ?n,j*z)|  <  TT  n  +  »;n'^»llzll)  <  ej=1  <  eM"Z*  . 

n  3-1  j-1 

Since  we  know  that  fn  is  uniformly  bounded  on  compact  subsets  of  £s  and 

f  (Z)  *  f(z)  pointwise  on  {zy  :  1 z |  <1,  z  e  t,  y  e  w}  the  convergence  is  uniform 
n 

on  all  compact  subsets  of  |fs.  This  completes  the  proof. 

An  obvious  corollary  of  this  result  is  the  following  characterization  of  the  class 
ES(W). 

Corollary  1.  An  entire  function  t  _i_s  W  -  lineal  if  and  only  if  it  is  in  Eg(w). 

An  essential  point  of  view  in  our  previous  presentation  was  the  analysis  of  an 
entire  function  on  0s  along  an  arbitrary  ray.  The  form  of  the  resulting  univariate 
function  allowed  us  to  make  conclusions  about  the  multivariate  entire  function.  On 
this  basis,  one  might  be  tempted  to  conclude  that  a  function  f  is  in  Eg  if  and  only 
if  f(zy)  is  in  Ei  for  all  y  e  Fs.  This  is  unfortunately  false  as  can  easily  be 
seen  by  the  following  simple  example. 


-7- 
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Example  2_.  Let  =  z*j  7^  -  1.  Then  along  every  ray  (ty^ty^),  f  is  the 

quadratic  polynomial 

Hty^ty.,)  =  t2(y^  +  y2?)  -  1 

i ch  vas  two  real  zeros  a*  1  so  is  in  F^.  However,  f  is  not  in  F2  because  its 
rer.i  set  is  the  unit  circle  (elements  ir.  E2  vanish  only  on  lines). 

This  example  can  be  substantially  generalized.  Let  q(-z)  be  any  (univariate) 
one-sided  Polya  frequency  function  given  by 


q(  z)  -  e  YZ  JI  (1  -  a .  z) ,  7,a.  >  0,  £  a.<°° 

1  3  ^  j=i  ~ 


and  define 

(2) 


f(Z)  =  q(z^+* • *+z2 ) ,  Z  *  (z  ,...,Z  )  . 

is  Is 


Then  it  is  cleaj  that  ^f(zy)  is  in  for  all  y.  However,  the  only  radial  funetio 

-z  -•••-z 
1  s 

in  Es  is  e  while  f  defined  by  (2)  is  always  radial.  We  will  come  back 

to  this  class  of  functions  in  the  last  section  of  the  paper. 

Let  us  also  point  out  that  there  are  functions  other  than  those  given  by  (2)  whic 

are  in  along  every  ray.  For  instance, 

2  2 

f  <  Z)  ~  z,  +  •  •  •  +  z  i  -  z 
1  s-1  s 

has  this  property.  The  analytic  form  of  all  such  functions  is  unknown.  However,  when 
rays  are  replaced  by  lines  the  problem  has  a  satisfactory  answer. 

Corol  lary  2_.  funct  ion  f  def  ined  on  is_  _i_n  Eg  if_  and  only  if  f  ( x+zv )  _is  J_n 

K 1  for  all  x ,y  e  Fs. 

proof .  If  f  f  Es  then  it  easily  follows  fro"  'definition  1  that  f(x+zv)  t  F ^ .  ?e 
prove  the  converse,  let  us  first  n<nte  that  since  for  all  x,v  t  F* ,  f(x*zy>  r  F 1  it 
follows  that  f  is  a  real  analytic  function  on  rC;. 

Tn  the  case  that  f(x)  *  h  t  x  f  P*5 ,  then 


f  l  x+?y ) 


-az 


for  some  real  constants  a  =  a(x,y)  >  0,  8  =  B(x,y).  Since  f  is  entire  we  may 


-AZ*Z+b*Z 


differentiate  -is  equation  to  see  that 

f(z)  =  e 

where  A  an  s  x  s  positive  semi-definite  matrix  and  beRs.  In  the  general  case,  we 
know  according  to  (3)  that  the  zero  set  of  f  consists  of  real  hyperplanes  r?  *z  »  t.. 
where  ^  e  RS,  lie  ■'ll  =  1  and  tj  £  Rs.  Thus  we  may  express  f  as 


,  m  n  j 

f  ( z )  =  e-A2"2+b-z  |  |  (z3.z>  lim  |  |  a  -  l)  , 

j=1  n+“  j=1  j  3 


where  E(z;p)  are  the  Weierstrass  polynomials  given  by 


E(  z;  0 )  =  1  -  z 


E  ( z ;  p ) 


(1  -  z)e 


P  =  1,2.. 


and  a^  are  any  integers  such  that  Z  ('j |  )  ”  <  “  for  finite  1 ,  [3,7].  There¬ 
'S  n 

fore  f(zy)  vanishes  only  at  z  =  5  ’y/t^  and  since  ftzy)  e  E,  it  follows  that 

i  (4^2  <  -  • 

j 


■  2 

Thus,  Zt  <  "  and  so  we  may  choose  a  =2.  Consequently,  f  e  E_  which  completes 
1  n  s 

the  proof  of  the  theorem. 

There  is  a  similar  theorem  for  the  class  Es(w)  where  W  is  an  admissible  wedne 
which  we  state  below. 

Corollary  3_.  A  function  f  f  Es(W)  _i_f  and  only  if 
1)  f ( x+zy)  e  E-,  for  all  x,y  e  Rs 

2  )  f  (  zy )  _is  _a  one-si  ded  Polya  frequency  function  for  al  1  y  f  w . 

PROOF.  Clearly,  accordinq  to  Definition  2,  properties  1)  and  ?)  hold  whenever 
f  e  ES(W).  Conversely,  according  to  the  Corollary  2,  property  II  above  implies  that 


-g- 
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Section  3.  LIMITS  OF  M7LT IV AFT ATF  P-*rrLT?:FR 

In  this  section  we  turn  to  tte  principal  apnl  i  c-atior.  r.f  the  previous  results  on 
entire  functions.  We  bee;;  p.  by  recaltirg  the  definition  of  the  multivariate  P-splin«. 
Definition  4^.  Let  , . .  •  ,xn  c  Fs  and  suppose  (x^, . .  .  ,xnl  =  convex  hull  of 

x°,...,xn.  T_f  fxC,#.../xn]  has  dimension  s  then  the  linear  functional 

*  '  0  n 

f  -  f(v  x  +  •  •  •  +  v  x  )  dv„  •  •  • 

'  0  n  1  n 


where 


S  -  {  (  V-  /  . . »  ,  V  )  :  V  *  n  ,  V  .  =  1  } 

0  r.  ■»  “  1 

0 


has  for  n  >  f  a_  (  un  in' jo)  contj  nuous  densi  ty  M(  x  !  ,  , .  .  ,  xr  ) 

n 


)  dv 


•'f  xl  V  )f(xHx  , 


b  great  deal  is  v.r.ow  .iv»it  t'-ic  f  un-t  \  or.  an  i  eert.j  :  n  linear  combination?*  ■  r  «•  i 
functions .  t suf  r  i  rw’,<‘  for  !>■;*  ;•»*••«  i  i  r- i  nr*  hr*r^  t*-at  w^pn  «  ••*  1 

t  i  on  co  rr  ecr-' *°  *,f,l  )  --c.-avn  v  jat  o  n-t.'-i  )*'•■.  it  j  <-  ..  •;(' 


0 


co  r  r 


polynomial  of  degree  n  -  1  with  n  -  2  continuous  derivatives  when  x®,...,xn  are 
distinct  (these  are  real  numbers  in  this  case).  In  general,  it  has  been  shown  that 
M(x|xn , . . . ,xn)  is  a  piecewise  polynomial  of  total  degree  <  k  with  k  -  1  continuous 
derivatives  when  the  convex  hull  of  every  subset  of  s  +  1  points  of  {x^,...,xn} 
forms  an  s-dimensional  simplex. 

In  this  section,  we  will  identify  all  distribution  functions  which  are  limits  of 
distributions  of  the  form 

X  X 

„  ,  ,  .  r  1  f  s  ,  .  0  ,  n  n ,  n .  , 

Fn(x>  =  n!  J  •••  j  M(y)x  ,...,x  ) dy 

.00  .00 

x  =  (x1,...,xs)  (we  will  sometimes  refer  to  Fn  as  a  B-spline  distribution).  Our 
basic  tool  in  what  follows  is  the  Radon  transform  of  a  distribution.  We  review  below 
the  elementary  facts  which  we  require. 

Let  u  stand  for  a  probability  measure  on  Rs.  Sometimes  we  will  deal 
simultaneously  with  probability  measures  on  R1  and  Rs.  However,  the  context  should 
make  the  domain  of  definition  of  the  measures  apparent. 

Definition  5_.  The  Radon  transform  of  ^a  probability  distribution  y  on  Fs  i_s  map 
from  (Is  =  {A  j  A  e  RS,  II A II  =  1}  into  probability  distributions  on  R^  defined  by 

(Ry)  (I  ,A)  =  u({y:y  e  RS,  A»y  e  I }  ) 

for  any  interval  ICR1. 

Sometimes  we  will  denote  the  Radon  transform  of  u  by  u ( * , X ) . 

1  _  s 

Lemna  _1/  For  all  t  e  F  ,  A  e 

f  “ito.  .  .  .  r  -it A  *y  . 

J  e  du(a,A)  =  J  e  du(y)  . 

1  s 

P  P 

rn 

PPOOF .  Let  g(t)  =  i  ft)  any  simple  function  on  p\  with 

i=1  1  i 

I.  "'I,  -  $ ,  j  *  Then 
j  l 


b)  If  y^t'.X)  +  y(*,X)  then  clearly  ( y )  +  T(y)  where  ?(y)  =  <J<  ( t ,  X ) ,  the 

characteristic  of  y(",X)  and  y  =  tX,  X  e  12s.  Our  hypothesis  implies  that  f(y)  is 

continuous  which  insures  that  it  too  is  the  characteristic  function  of  some  probability 

measure  y  on  Rs.  (c.f.  H.  R.  Pitt,  Integration,  Measure  and  Probability,  Oliver  and 

Boyd,  London  1963).  Thus  by  Lemma  1  (Ry)(»,X)  =  y(*,X)  and  y  *  u. 

n 

We  will  say  that  a  probability  measure  y  on  Rs  is  a  Polya  distribution  if  the 
reciprocal  of  its  Laplace  transform  is  in  the  (normalized)  Polya-Laguerre  class  on 
Rs,  that  is, 


/  e  Z  Xdu(x) 


_ 1_ 

f(z> 


where 


f(z)  -  e**2*2*5 


’Z  I 

3=1 


(1  +  C^‘Z)e 


-^•Z 


(It  should  be  understood,  if  not  otherwise  explicitly  stated,  that  the  range  of 

validity  of  the  above  expression  is  the  hyperstrip  Re(C^’Z)  >  -1  in  £s.) 

A  B-spline  distribution  y  has  the  form 

y  (I)  =  ni  /  M(x|x°'n, . . . ,xn,n)dx 
n 


for  some  x®'n 
Proposition  2_. 
distributions . 
PROOF.  Since 
(3) 


. . . ,xn'n  f  Rs. 

If  y  _is  a^  Polya  distribution  then  it  is  a_  limit  of  P-spline 


y  is  a  Polya  distribution  we  have 

/  e-Z,ydy(y)  = 
RS 


where  f  is  a  normalized  Polya  freauency  function.  Thus  accordina  to  Motzkin  and 
Schoenberg  (see  the  remark  following  Theorem  1)  there  is  a  set  of  vectors 
(x°'n . xn'n)  in  Ps  such  that 
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1  i  pi  j  I  +  z  •  x J ‘  ‘ )  -  *  <  z  ) 

n-f-t>  -  _r  n 


uniformly  on  compact  subsets  of  Re«  Thus  for  everv  X  f  I'-'" 


f  (  :•  A )  --  Mb  I  !  (1  *■  7\-x~'n) 

r-"Mr  i~C< 

uniformly  on  compact:  subsets  of  R^.  Hence  by  the  Curry -Schoenberg  theorem  we  know 
that  the  secjuer.ee  of  distributions 

f  4 )  u  (I,).)  =  n!  ;  J«(jlX‘x,'n,...,X»xr*'n)rto.  I  C  P1 

*  I 

converges  to  a  distribution  u(t,\)  which  satisfies 


(5) 


e"z^,j(o/A) 

1 


1 


f  (?.\1 


According  to  (3),  (Pu)(*#X)  -  Moreover,  if  we  define 

(6)  -u  <I>  =  ^xlx°'n . xr'n)dx,  I  C  ps 


then  we  can  see  that  the  Falon  transform  of  i.  (I)  is  u  (I.X ),  because  by  Definition 

o  r. 

4  and  Lemma  J,  if  I  P1 

(Pi.  )(I,X)  =  f  (  *  •  X )  V  ^  x  I  X  '  n,  .  .  .  ,x  1  )  dx  *  .  y,  <Zv  .A*x^  ,ri)  dv  •  •  *nv 

n  I  J  AT  ■>  In 

P  c’! 

.  f  r*i  r-  *-> 

-  r. !  v r  t !  W  X  -v  ' '  )  ■»:.  --  u  ( I ,  X  )  . 

:  n 

Thus  we  have  sbovn  that  (K.  M*,\)  (  rw)  .  Accordina  to  Pronosition  1,  y  ♦  y 

n  ‘  n 

because  the  characteristic  ‘unction  of  i  s  obviously  uniformly  continuous  in 

a  neighborhood  ->f  re.  This  lutes  ?  w.  -rroof . 
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Proposition  3_.  Suppose  y  jus  a^  probability  measure  on  Rs  and 

y  (I)  =  n!  /  M( x |x° ' n1 . . . ,xn,n) dx 

n  i 

converges  to  y.  Then  u  _is^  a_  Polya  distribution  with  Laplace  transform 

j  e  *dy(x)  =  —  -  - 
R® 


where 


n  j  ,n# 

lim  T“T  (l  +  - - — )  -  f  (  z)  . 

n~  5-1 


PROOF.  As  before  we  have 

(Ry^HI.X)  -  n  1  /  M(t  |  X*x° .  ,X*xn  ,n)dt 

and  since  (Ru  )(X»X)  *  (Ry)(I,X)  for  each  X,  the  Curry-Schoenberg  theorem  implies 
that 


n  j,n 

lim  |  |  (l  +  z  — - )  -  g(z,X) 

n-*»  j-0  n 


converges  uniformly  on  compact  sets  for  every  fixed  X  to  some  g(*,X)  e  e  •  Mow  just 
as  in  the  proof  of  Theorem  1  we  conclude  that 

g(z,X)  =  f(zX) 


where  f  is  given  by 


n  j  ,  n  m 

(7)  lim  TT  (l  +  * - — )  =  fC)  . 

n-*®0  j=0  n 

This  convergence  is  uniform  on  any  compact  nart  of  any  ray  through  the  origin.  The 
Curry-Schoenberg  result  tells  us  that 


-IS- 


/  e  ita<3;j(o,A) 

p1 


'  1 

gTtTXT  “  TTtTT 


and  so  invoking  Lemma  1  again  we  obtain  the  desired  result,  namely 

/  e"'’Ydufy>  =  — . 

RS 

To  complete  the  proof  we  only  have  to  show  that  the  convenience  in  (7)  is  uni  for" 
compact  sets  of  (fs. 

In  (7)  the  derivatives  of  the  product  converoe  to  the  respective  derivatives 
f.  Thus  we  see,  as  in  Theorem  1,  that  these  products  are  uniformly  bounded  on  an 
compact  set  and  so  the  convergence  is  uniform. 

Propositions  2  and  3  combine  to  give  the  main  theorem  of  this  section. 
Theorem  _3 .  _A  probability  distribution  u  j_s  the  limit  of  B-spline  distributions 
and  only  if 

,  e  ydu(  y)  »  ypj 
Ps 


where  f(0)  =  1  and  f  e  Eg. 

Corollary  4_.  necessary  and  sufficient  condition  for 

U  (I)  =  n!  /  M(x|x°'n . xn'n)dx 

n  * 


to  converge  to  the  distribution  y  defined  hy 


dii(y ) 


i s  that 

a)  1  im  —  max  !. x  ' ,n  i:  -  h 

n-»“>  najyn 


-If- 


lim  -  l  xi,n  =  y 

n-»«  n  j=o 


Example  £.  Let  xr  m,  0  <  r  <  m,  be  any  scalar  sequence  such  that  the  corresponding 

-1  /2  -t2 

univariate  B-spline  distribution  converges  to  the  Gaussian  distribution  Tt  e  dt. 
(The  conditions  on  xr  m  needed  to  insure  this  is  given  by  Corollary  4.)  Now,  form  a 
vector  sequence  x^'n,  0  4  j  <  n  =  ms,  by  "placing”  xr  m  on  each  coordinate 
direction.  Then  we  may  again  use  Corollary  4  to  conclude  that  the  corresponding 

multivariate  B-spline  converges  to  the  multivariate  Gaussian  distribution 

-8/2  -«xll2 
it  e 

We  will  say  that  a  B-spline  distribution 

y  (I)  =  n!  '  M( x |x° 'n, . . . ,xn ,n) dx 
"  I 


is  W-restricted  if  x^'n  e  w®. 
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Theorem  4^.  prohabi 1  ity  distribution  u  is  the  limit  of  W-restricted  B-spline 

distributions  p  if  and  only  if 
— - - n - • - 

(8)  /  e-’-Vauty)  =  TiTT 

s 

R 

where  f(0)  =  1  ana  f  e  E  (W). 

We  will  prove  this  theorem  in  two  propositions. 

Proposition  4..  u  xs_  a_  distribution  satisfying  ( 8 )  then  It  is  a_  limit  of  w- 

restricted  B-spline  distributions  and  supp(u)  C  cohull (W°). 

PROOF.  According  to  Corollary  1,  there  is  a  sequence  {x^ ,n, . . . ,xn,n}  C  W®  such  that 


lim  |  |  (1  +  z*x^'n)  =  f(z)  . 
n+“  j=0 

Thus,  as  in  Proposition  2,  u  (I)  --  n!  /  M(  x |x° ,n, . . . ,xn ,n) dx  converges  to  u,  since 

"  X0  n  n  n 

f  e  E  .  Moreover,  because  supp  u  C  [x  '  ]  it  follows  that  the  support  of 

s  n  — 

U  C  cohull {vP  ) . 

Proposition  5_.  Suppose  u  _is  a^  probability  measure  on  Rs  and 

Un<I)  =  n!  /  M(x|x°'n . xn'n)dx  , 


where  xJ'n  e  W*"1.  If  u  conv»roes  to  u  then 

— “  -  n  ..  r  n-r  — —  "r_ri  _T 

distribution  with  Laplace  transform 

r  *-  Z  »y 
j  e  ydp(y)  * 


U  j_s  a_  W-restricted  Polya 


_1 _ 

f(z)  ' 


n 

f(z)  =  lim  [  f 

n+co  j=n 


an  1  supp(*j)  2.  cohulW  )• 

This  result  is  proved  hy  o'"iu"'°nts  similar  to  those  used  to  prove  Proposition  4.  We 
omit  th*»  details. 


Cl  _ 


1 


The  analog  of  Corollary  4  is 

Corollary  5_.  Let  W  be  _an  admissible  cone .  Then  a_  necessary  and  sufficient  condition 
for 


U  (I)  *  nl 
n 


■  M( x | x^ ' n , . . . ,xn ,n ) dx, 


W 


0 


0  0 

to  converge  to  the  Dirac  measure  at  5  e  W  _is_  that 


a) 


lim  —  max  llx^'nll  =  0  , 
n*"  n  0<j<n 


b) 


..  1  V  j,n  a 

lim  -  >  xJ  -5 

n+“  j=0 


Section  4.  INVERSION  FORMULA  AND  SMOOTHNESS  PROPERTIES  OF  POLYA  DISTRIBUTIONS 

In  this  section  we  investigate  some  further  properties  of  Polya  distributions. 
Their  characteristic  property,  as  defined  earlier,  is  that 


/  e-2*ydu( y)  -  Re(;j*z>  >  -l 

FS 


where 


f(z) 


-AZ*Z+C 


I  I  (1  +  «Z)e 

j=1 


£  (J^ll  <  “,  A  is  a  positive  semi-definite  s  x  s  matrix,  and  r ^  e  ps, 

j=1 

j  3  • 

We  begin  with  the  finite  kernel  case,  i.e.. 


(9) 


/  e"^*yau(y)  =  — - - -  ,  Pe(Cj*z)  >  -1,  e  Rs  . 

PS  TT  (1  + 


-•n- 


In  general,  we  can  always  represent  this  measure  as  an  n-dimensional  integral  in 
the  following  way.  We  define  the  linear  functional 


»  oo  -( t  +  • 

T  (4>)  =  /•••/  e  1 

o  o 


•+t  ) 

n  4>(t,;  +...+t  ;n)dt 
i  n  i 


•dt 

n 


and  observe  that 


T 

C 


1 


(e 


-(v 


•+t  ) 
n 


V 


•dt 


1 


TT  (i  ♦  2*53) 

j-i 


where  Re(z*C3)  >  -1,  j  «  1,...,n. 

Thus,  for  u  defined  by  (9),  we  obtain 

T(f)  -  /  f  ( x) du( x) 

»S 

at  least  for  functions  which  are  bounded  and  continuous  on  Rs.  Perhaps,  it  should  be 
pointed  out  that  this  "Inversion  formula"  for  the  measure  (j  in  the  univariate  case  is 
not  the  representation  given  by  Hirshman-Widder ,  [4].  A  “conventional"  representation 
would  express  U  as  an  s-dimensional  integral.  To  this  end,  let  us  observe  that  if 
m  »  dim  span(c1 , . . . ,£n}  <  s  then  T  does  not  come  from  an  b’-density,  i.e.  there  is 
no  L’-density  G(x)  such  that 

T(f)  «  /  f (x)G(x)dx  . 


Otherwise,  there  would  be  a  X  e  Fs  -  {0}  such  that  X*;-'  =  0,  j  =  1,...,n.  Then 


.  m,  ”itX»x,  r  -itX*x 
1  =  T(  e  )  =  j  e 


i 


G(x)dx  . 


r 


However,  because  G  e  lVr®)  the  limit  of  the  right  side  of  this  equation  is  zero  for 

t  +  ®, 

When  m  •  s  we  consider  two  cases.  For  n  «  s  we  write  «  Ae-’ ,  j  =  1,...,n 
where  (e^>k  “  3'k  *  1,..««n  and  A  is  the  nonsingular  s  x  s  matrix  whose  j-th 

columns  is  the  vector  Thus,  in  this  case, 


1  -1MA_1)TX  1  B, 

Idet  A|  ,  x  e  [(  . C  Jc 


(10) 


G(x)  » 


otherwise 


Is  Is 

where  [(,..»,(]  is  the  cone  spanned  by  5  ,...,;  • 
c 

The  last  and  principal  case  occurs  when  tn  =  s  <  n.  In  this  case  we  define  the 

density  corresponding  to  T  recursively.  We  suppose  for  ease  of  notation  that 
1  8 

C  ,...,?  are  linearly  independent  and  define 


(ID 


G(xU 


/  e_tG(x  -  tCi+1|C1,... ,C*)dt 
0 


initialized  by 
(12) 

where  G(x)  is  given  by  (10).  It  is  an  easy  matter  to  prove  that 


GUI?1,. ..,5s)  “  GU> 


(13) 


S 


(f)  -  /  f(x)G(x!s\...  ,Sn)dx 


(induction  will  do).  Since  T  .  is  invariant  under  a  permutation  of  the 

,  ’  r" 

S  ,•••,<» 

1 

vectors  s  . 


Sn  so  too  is  Gixls1 , • • • ,?n) •  Consequently,  the  recursion  above  can 


begin  at  any  set  of  linearly  independent  vectors  and  proceed  throuah  any  ordering  of 
the  remaining  elements  of  k’ , . . . .  We  might  also  mention  that  this  fart  can  be 
used  to  bound  G( x| C 1 , . . . ,;n) .  In  particular. 
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L 


m 


maxlG( x | £ 1 |  <  min  vol  'k 
x  K  S 

i1  is 

where  the  minimum  is  taken  over  all  s-dimensional  simplicies  K  =  [0,5  ,...,5  1 

formed  by  any  s  vectors  from  (c1,...,;55}  and  the  origin  in  P.s. 

Because  of  (11)  it  is  to  be  expected  that  G ( x| 51 , . . . ,5°  )  has  a  number  of 

continuous  derivatives.  To  make  this  precise,  we  denote  the  directional  derivative 

of  f  in  the  direction  y  e  Rs  by  D^f,  Then  (11)  clearly  implies  that 

(14)  (D  +  DG'xIc1 . 5i+1)  =  G(xlc’,...  ,5*> 

^  IS. 

at  every  point  of  continuity  of  G(x|5  ,...,5  ).  In  general,  we  see  that 


(15) 


D  G  ( x  |  5 

y 


.,cn) 


l  u.Gtxk1 - ,5j-V+1 . Cn> 

j-1  3 


where 


y 


j=i 


V 


i  v  -  o  . 

3»i  3 


The  proof  of  (15)  uses  (14)  and  the  fact  that  G ( x I  5 1 , . . . ,5n)  is  invariant  under  any 
permutation  of  the  vectors  51,...,5n.  Returning  to  (11)  we  note  that  G( x | 5 1 , . . . ,5S+1 ) 
is  continuous  on  Rs,  if  the  vectors  51,..>»5*  1 ,5*+ 1 » • • • ,CS ,  5S+ 1 ,  i  =  1,...,s,  are 

Q  54"  1 

linearly  independent.  Since  in  this  case  for  every  x  e  Rs  the  line  x  -  t5 

1  *5 

intersects  the  cone  [5  ,...,5)  at  most  once.  Thus  we  may  combine  our  previous 

c 

remarks  to  conclude  that  whenever  every  subset  of  s  vectors  from  {5 *,...,  5n }  is 

linearly  independent  G  (  x  I  5 1 , . . . ,  rj' )  has  n  -  s  -  1  continuous  derivatives. 

In  one  dimension  the  precise  form  of  G  ( x!  51  , . . . ,  5n )  as  a  combination  of  two 

exponential  polynomials  is  known,  Hirshman  and  widder  (Theorem  8.2,  page  31)  [4],  For 

the  ourpose  of  presenting  a  similar  result  in  Ps  we  prove  the  following  lemma. 

1  n 

The  points  C  /•••,;  are  saul  to  be  in  general  position  if  every  subset  of 


1  noint*s  forn  r.^nni  1  si^pl^x 


Let  be  in  general  position.  Then  for  every  set  I  C  {l,...,n} 

containing  III  =  s  elements  there  is  a  unique  Xj  e  Rs  such  that  1  +  =  0. 

j  e  I.  Moreover,  Xj  necessarily  has  the  property  that  1  +  £3*xx  *  0,  j  i  I. 
Lemma  2 .  Let  {0 , 5  \  . . .  , i;n}  be  in  general  position.  Then  for  every  polynomial  Q 
total  degree  <  n  -  s  we  have 


Q(x) 


I  Q<x  ) 

IH=s  1 


I  !  (1  +  ;j*x) 

in _ # 

TT  (i  +  cj*x  ) 

in 


PROOF ■  The  space  of  polynomials  of  total  degree  <  n  -  s  has  dimension  N  *  (n). 
Thus  it  suffices  to  observe  that  the  N  polynomials  (of  degree  <  n  -  s) 

TT  (1  +  S3  *x) 

Q  (x)  =  ^ - - - ,  IC  {l,  . . .  ,n}  ,  III  =  s 

I  |  (1  +  C3  *x  ) 
j/I 

satisfy  Qj(Xj)  =0,  X  *  J  for  then  they  must  be  linearly  independent. 

We  may  now  specialize  Lemma  2  to  the  polynomial  Q(x)  =  1  to  obtain  the  partial 
fraction  decomposition 

1 „  v  a  _ 1 

TT  (i  +  IT=s  1  TT  (i  +  c3-z) 


where 


I  I  (1  +  C j*x  ) 

in 


Hence  we  obtain  the  basic 
(16)  Glxlc' 


relation 


il|=s 


axG(x!c 


I  -  U. 


i 

s 


}  . 


of 
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/  V*!*:1 . cn)dx 


*»  ,  u 

*  lim  /  f(  [  v.C3}uj(  y  v.)dv  •••dv 

n+O  n  j*=1  3  j»1  3  n 

y  v  ,<h 

j-i  3 


«  lim  hn  /  f(h  ^  y  v.C^jufh  ^  y  v.)dv  •••<5v 

h*0  n  j»1  3  j=1  3  1  n 

y  v  <i 

j-i  3 


h+0  0  n 

y  v  -T 

j-i  3 


lim  h  n  /  (  /  f(h_1  4  v.?3)u(h  'ijdv  **»dv  ;dT 

v.^A  n  _  .i_«  3  ■  n”  ’ 


lim+  h  n  /  Tn  u>(h  t)(/  f(h  T  £  v . )dv1  •  •  ) dr 


h+0 


,n*1 


j-1 


u  *  .  . 

“  lim+  /  /  <u(h)hn  s  f(x)M(h  x|C  ,...,Cn)dx  . 

h*0  Rs  0 


Thus  we  obtain 


G  (xk1 . ;n)  -  /  aXh)hn'S'1M(T'1x|e1 . Cn)dx  . 

“  0 


Besides  the  case  u(h) 


which  gives  the  formula 


„ ,  ,1  n,  r  -n.n-s-i  ,  -1  _n,  . 

G(xC  ,...,5  )  =  J  eh  M(T  xU  «•••,*,  )  dt 

0 


the  choices  u(t)  *  Xrn  and  w(t>  =  1  yield  the  known  formulas 

l  0  /  ■  J 


M(x|o,c1,...,;n)  =  /  t  n+s  %( tx| .  ,?n)dT 


and 
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h(x | ? 1 , . . .  ,;n)  —  /  t  n+p  Vtixlc1 . ;n)d-c  , 

o 

respectively,  [6!.  (This  last  identity  is  true  when  0  t  U 1 ,  . . . ,£n] . ) 

Next#  we  turn  to  the  infinite  kernel  case.  For  that  purpose  we  introduce 
Definition  7_,  function  f  _in  Eg  given  by 


0  «  j 

f(  r)  =  e  ‘11(1+  ^  ' 

j-1 


is  called  degenerate  if  and  only  if 

S  —  span  t a  ^ , . « . , a  , £  < r,  ^  K 

where  a^  are  the  rows  of  the  matrix  A. 

We  will  also  say  A ( x ) ,  x  e  Ps ,  i_s  ^  Polya  frequency  function  provided  that 
du( x)  =  A(x)dx  is  a  Polya  distribution. 

Theorem  5,  Let  f  e  Eg  and  suppose  f(0)  =  1.  Then  1/f  is  the  Laplace  transform  of 

A  po*Va  frequency  function  if  and  only  if  f  is_  nondegenerate. 

PROOF.  Suppose  f  Eg ,  f ( 0 )  -  1 ,  is  degenerate  and 

=  /  e~‘*XA(x)dx,  Re(Cj*z)  >  -1,  j  =  1,2,... 

'  Rs 

for  some  density  function  A(x). 

Since  f  is  degenerate  there  is  a  A  e  Rs ,  A  *  0,  and  A  i  S.  Hence 

zr®  .  A 

f(Az)  ■  e  and  so 


0 

e'lt? 


•  X 


f  -itX  #x. 
j  e  A(x)dx 

s 

P 


for  all  t  e  P1 .  This  is  obviously  a  contradiction  because  the  right  hand  side  of  the 
equation  goes  to  zero  as  t  ♦  00 ,  by  the  Pieman n-Lebes cue  theorem* 

Conversely,  suppose  f  is  nondegenera .  We  will  now  show  how  to  construct  a 
density  whose  Laplace  transform  is  1/f.  First  let  us  observe  that  if  f 1 ,  f2  are 


-26- 


I 


in  Eg  and  1/f1  is  the  Laplace  transform  of  a  density  A(x).  Then  1/f,  for 

f  =  fif2'  is  the  LaP^ace  transform  of  the  density 

/  A(x  -  y)du(y) 

Rs 

where  du  is  the  Polya  distribution  corresponding  to  fj.  Thus  to  prove  the  theorem 

it  suffices  to  find  a  factor  of  f  which  corresponds  to  a  Polya  frequency  function. 

When  A  is  of  full  rank,  r  =  s  then  the  well-known  formula 

eAz*Z  .  - 1 -  /  e-Z*VU"1x'x)dx 

(2/if)Si/det  A  rs 

already  provides  the  density.  In  the  other  extreme  case,  when  A  »  0,  there  are  s- 

linearly  Independent  vectors  from  :  j  =  1,...}  say  t  ,  ...,t  .  Then  the  formula 

1  s 

1  r  —  2*  x 

- -  “  /  e  G(xlc1  ,...,i;s)dx 

TT  (1  +  ?J*2>  RS 

j-1 


which  we  used  earlier  in  the  section  gives  us  a  density.  Finally  if  0  <  r  <  s  then 

1  £  1  Jt 

there  are  l  -  s  -  r  vectors  say  C  ,...,£  such  that  {a1,...,af,(  ,...,C  )  spans 
Fs.  In  this  case,  the  function 

l 

F(z)  =  e_AZ'Z  11(1+  Cj’Z> 


is  a  factor  of  f.  To  see  that  it  corresponds  to  a  density  we  may  as  well  make  a 
linear  change  of  variables  and  assume  it  has  the  form 


2 

•  «-z 


r 


I  I 

J-r+1 


(1 


V 


Now,  it  is  clear  that 


_ 1_ 

F  (  Z ) 


- - —  /  e-Z*XA(  x)  dx 

(2/i)z  Rs 
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where  A  is  the  density 


This  last  case  proves  the  theorem. 

Let  us  end  this  section  by  noting  that  the  density  A(x)  will  be  C^tR3) 
provided  that  either  A  is  positive  definite  or  that  there  is  a  set  (c3:  j  e  j) , 
|J|  *  which  are  in  general  position  with  the  origin  in  Rs.  This  latter  fact 


follows  from  our  earlier  discussion  of  the  finite  kernel  case. 


Section  5.  POLYA  FREQUENCY  FUNCTIONS  AND  VARIATION  DIMINISHING  TRANSFORMATION  ON  RS 
In  this  section  we  will  indicate  the  extent  to  which  we  can  generalize  the 
variation  diminishing  property  of  univariate  Polya  frequency  functions  to  higher 
dimension.  Our  remarks  are  based  on  the  following  elementary  observation. 

A  ridge  function  has  the  special  form 

f^(x)  «  g(w*x) 

where  g  is  a  univariate  function  and  »  e  0s,  the  unit  sphere  in  Rs.  The 
convolution  operator 

(A*f)(x)  »  /  A(x  -  y ) f ( y ) dy 
RS 

maps  the  ridge  function  f^  into  the  ridge  function  h(u*x)  where 

h<t)  •  /  A(y)g(t  -  u»y)dy  . 


The  induced  map  T  defined  by  h 
w 


where 


(R  A) (t)  -  /  A ( x ) dx 

X*U“t 


is  the  Radon  transform  of  A 


Theorem  5^  _lf  A  _is  _a  Polya  Frequency  function  on  Rs  then  for  any  a  e  f!8  the  man 
Tu  defined  a*,  ve  Is  variation  diminishing. 

PROOF.  To  prove  this  re.ult  we  must  demonstrate  that  (R  A)(t)  is  a  Polya  frequency 


function  in  t  for  every  u  f  !1  .  Since 

r  e_2,y 


A(y)dy  = 


f  (  Z ) 


where  f  e  Eg  we  see  that  the  Laplace  transform  of  R^A  is  f  ^(zw«y).  But  f(za»y) 
is  in  the  class  •  This  implies  by  one  of  Schoenberg's  basic  results  17,8]  that 
(R  A) ( t)  is  a  Polya  frequency  function. 

The  above  result  only  gives  a  sufficient  condition  for  A  to  be  a  Polya  frequency 
function  on  Rs.  There  is  an  extensive  collection  of  radical  functions  satisfying  the 
conclusion  of  Theorem  5  but  which  are  not  Polya  frequency  functions.  Specifically,  we 
let 

A(x)  »  h(Hx»2),  II x II 2  »  x2+»»»+x2, 

where 


„  00  J,  ,„(to) 

1^#*.%  1  f  S/2  I  S-2  )  / 2 

^  s/2  (s-2)/2  '  0  2  ' 

(2t)S/2t(  0  q(-J  > 


(the  Bessel  transform  of  l/g(-o  ))  and  g  has  the  form 


g(t)  =  e’Yt  II  (1  +  at), 
j-1  3 


y  >  n. 


j-1 


a  <  00  . 
1 


It  is  well-known  that  the  Fourier  transform  of  A 

/  e  iX  yA(y)dy  = 


is  given  by 
a  (  -  n  x  I! 2  ) 


and  thus 
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/  e  Z*yA(y)dy  - - - - - - — ,  z«  (e1>...,zfl)  . 

R8  9*21  +  *  ’  "+zs  ) 


Now,  it  quickly  follows  as  in  the  proof  of  Theorem  S  that  (R  A)(t)  is  a  Polya 
frequency  function  in  t  for  every  ui  e  fl8. 

Postscript 

In  their  paper  Curry  and  Schoenberg  showed  how  the  Brunn-Minkowski  Theorem 
concerning  the  volume  of  a  convex  body  intersected  with  a  hyperplane  leads  to  a 
"convexity"  inequality  for  the  univariate  B-spline.  The  same  argument  gives  similar 
information  about  the  multivariate  B-spline,  We  record  this  result  below. 

Recall  that  a  nonnegative  function  g  defined  on  P1  is  a  Polya  frequency 
function  of  order  two,  if 


»(W 

*(W 


9(W 

g(t2-o2) 


>  o 


whenever  t1  <  . 

Theorem  T_.  For  any  x,y  e  R8  the  function  g(  t)  ”  M(  x  +  ty|x°, . . .  ,xn)  _is  a_  Polya 
frequency  function  of  order  two. 

PROOF.  Using  the  geometric  interpretation  of  the  multivariate  B-spline,  [2,6],  as  the 
volume  of  a  polyhedra  and  the  Brunn-Minkowski  theorem  it  follows  just  as  in  [1]  for  the 
univariate  case  that  (M(x  +  ty|x°, . .. ,xn) ) '/n-s  is  a  concaVe  function.  The  theorem 
now  follows  by  observing  the  simple  fact  that  every  positive  power  of  a  nonnegative 
concave  function  is  a  Polya  frequency  function  of  order  two. 
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